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Summary 

The  goal  of  this  work  was  to  inveatigate  the  performance 
limitations  of  the  "standard"  GaAs  FET  structure.  During  the 
contract  period  we  have  constructed  an  ASCI3  epitaxial  system 
used  to  provide  buffer  layers  for  our  FET  structures,  we  have 
developed  a  submicron  lithographic  processes  using  deep  O.V. 
techniques  and,  employing  these  techniques  we  have  produced  work¬ 
ing  .5  micron  gate  devices.  In  addition,  we  have  investigated 
the  "gettering"  of  substrates  as  a  technique  to  improve  the 
mobility  of  ion  implantated  layers.  The  result  of  this  experi¬ 
ment  showed  was  a  correlation  between  improved  hall  mobilities 
and  gettered  substrates.  Finally,  a  investigation  of  a  self 
aligned  source  drain  structure  has  commenced  and  several  theoree- 
tical  studies  are  reported. 


Materials  Development 


In  section  A,  ve  will  report  on  the  progress  in  develop- 
ment  of  a  As C 1 3  system  for  Che  production  of  high  quality  buffer 
&  active  layers  for  our  microwave  devices.  It  is  anticipated 
that  the  ASCI3  material  will  provide  thick  buffer  layers  for  iso¬ 
lating  the  substrate  effects  while  growth  technologies  such  as 
MBE  will  provide  the  thin  critical  active  regions.  Part  B  of 
this  section  concentrates  on  our  gettering  experiments  on  GaAs 
substrates  for  direct  ion  implantation  applications  while,  part  C 
reports  on  the  status  and  capabilities  of  our  recently  acquired 


MBE  system. 


A 


Project  Stitut  of  CVD  Reactor 


Introduction 

tfe  report  here  the  status  of  our  redesigned  halide  trans¬ 
port  VPE  reactor  and  results  obtained  to  date.  This  effort 
parallels  construction  of  a  pyrolytic  CVD  reactor  now  completed. 

As  previously  mentioned,  GaAs  FET  performance  shows  a 
great  dependence  on  submieron  doping  thicking  products,  backround 
impurities  and  substrate  properties.  The  system  is  designed  to 
grow  thick  (greater  than  10  micron)  high-resistivity  buffer 
layers  which  isolates  our  FET  devices  from  substrate  related 
problems . 

The  achievement  of  high  resistivity  buffer  layers  has 
been  well  demonstrated  in  the  literature  to  require  use  of  the 
well  known  "mole  fraction  effect"  first  described  by  Dilorenxo 
and  Moore  in  1971.  A  two  bubbler  systems  similar  to  that  of  Cox 
and  Dilorenso  (1971)  is  adapted,  but  with  certain  changes  and 
simplifications. 


-  *  - 


VPB  Syatei  Technoloi 


A  systematic  outline  of  our  system  as  originally  desig¬ 
ned  is  shown  in  figure  (1).  It  was  at  first  intended  to  deaos- 
trate  high-purity  and  reasonable  process  control  with  a  simple 
reactor  design.  Systems  improvements  were  later  added  as  they 
proved  necessary. 

The  original  design  was  a  teflon  plumbing  system  using  a 
continuous  hydrogen  purge  to  avoid  adding  reactor  bypass  valves 
as  a  possible  source  of  contamination.  The  use  of  teflon  is  in¬ 
tended  to  minimize  the  presence  of  any  metallic  impurities,  as 
well  as  decrease  the  "memory  effect"  of  any  system  reagents. 

This  teflon  plumbing  and' buffer  valves  still  seem  to  be  a  useful 
system  feature  and  were  retained  in  use. 

The  first  growths  utilized  ASCI3  obtained  in  pre¬ 
packaged  plug-in  bubblers  from  .Apache  Chemicals,  Inc.  These 
ASCI3  bubblers  were  designed  to.  plug  into  solid-state  temperature 
controllers  also  manufactured  by  Apache  chemicals,  Inc.  This  was 
deemed  to  be  a  valuable  feature  which  would  minimize  toxic 
hazards  associated  with  loading  the  CVD  system  with  AsCl3,  ya 


•  ... 
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Flow  Schematic  of  GaAs  Vapor  Phase  Epitaxial  Reactor  System 


Figure 


experienced  many  hardware  problem*  with  these  temperature  con¬ 
trollers.  After  being  sent  several  updated  versions,  we  have 
only  this  year  obtained  models  which  work  as  prescribed. 

To  check  the  reactor  kinetics  45  grams  of  liquid  Gallium 
as  source  material  was  loaded.  We  then  attempted  source  satura¬ 
tion  at  820°C  and  observed  the  saturation  time  and  other  charac¬ 
teristics.  We  found  it  necessary  to  add  a  aource  baffle  to  in¬ 
crease  source  saturation  efficiency.  Growth  parameters  of  liquid 
source  runs  3  thru  6  are  shown  in  Table  I. 

The  source  boat  was  later  loaded  with  approximately  40 
grama  of  crushed  GaAs  as  a  solid  phase  aource.  Solid  sources 
have  been  found  to  offer  better  thickness  control,  without  the 
troublesome  source  saturation/etch  cycle  of  liquid  sources.  Bet¬ 
ter  surface  morphology  is  also  observed. 

Initially  we  used  as  source  material  the  upper  "cone” 
portion  of  Bridgeman  process  semi-insulating  GaAs  crystals  pro¬ 
vided  by  Cominco,  Inc.  This  solid  source  material  was  used  in 
growths  under  a  variety  of  reactor  temperatures,  carrier  gas 
flows,  and  bubbler  temperatures  shown,  in  Table  I,  runs  7  thru  21. 


Over  certain  parameter  ranges  good  morphology  and  thickness  con- 


V- 


crol  were  obtained.  Certain  growth  and  hardware  problems  were 
encountered  which  we  attacked  with  several  design  changes. 


r.v;.".' 
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We  decided  to  reduce  the  normal  buildup  of  reaction  by 


products  at  the  exit  end  of  the  reactor  tube  by  use  of  a  hydrogen 
counterflow  system.  This  build  up  could  produce  particle  con¬ 
tamination.  Several  endcap  designs  were  tried  in  which  the  sub¬ 
strate  holder,  dump  tube  and  endcap  were  fuzed  into  one  quartz 
containment  assembly  as  shown  in  figure  (2).  The  solid-source 
was  changed  to  sealed-vessel  synthesized  GaAs  dices  supplied  by 
Morgan  Semiconductor  for  the  runs  remaining  runs  after  run  #20. 
This  high  purity  source  material  gave  only  a  superficial  improve¬ 
ment  in  the  electrical  characteristics.  Subsequently  we  were 
able  to  eliminate  'the  solid  source  as  a  cause  of  the  extrinsic 
compensation  since  no  source,  temperature  or  other  growth 
parameter  dependence  was  observed. 

While  at  this  point  the  faulty  temperature  controllers 
were  still  being  used,  it  was  suspected  that  the  ASCI3  itself  was 
contaminated  and  causing  the  epilayers  to  become  heavily 
compensated . 


The  ASCI3  wa*  returned  in  the  plug-in  bubblers  to  the 
manufacturer  for  analysis.  Flame  spec trochemical  analysis 
revealed  approximately  50ppm  Fe ,  20ppm  Si  and  a  few  ppm  traces  of 
Mg  and  other  elements.  The  source  Fe  levels  are  very  consistent 
with  the  electrical  characteristics  observed  in  the  epilayers. 


Fe  backrounds  generally  in  Che  10"®  mole  fraccion  range  will 
yield  high  resiscivicy  layers  wich  mobilities  of  approximately 
1000-2000  cm^  V'^Sec-^.  We  tentatively  concluded,  therefore, 


that  our  low  mobilities  and  high  resistivities  were  caused  by 
unintentional  incorporation  of  bubbler  Fe  due  to  contaminated 
As  C 1 3  . 

System  Modifications 

To  alleviate  the  above  difficulties  we  have  made  changes 
in  our  source  materials  and  proccess.  The  most  important  of 
these  by  far  is  the  use  of  good  quality  ASCI3  obtained  from 
other  domestic  manufacturers.^  We  will  fill  our  plug-in  bubblers 
with  the  pure  As  Cl  3  demonstrated  in  the  field  as  having  the 
requisite  purity. ^  This  change  alone  is  expected  to  solve  most 
of  the  compensation  problem. 

The  reactor  tube  itself  has  been  modified  extensively. 

Ac  the  inlet  side  a  vacuum  chamber  thermally  isolates  the  inlet 
tubes  from  the  source  region.  This  aids  greatly  in  reducing  any 
premature  thermal  decomposition  of  main  or  bypass  injection  line 
reagents.  As  already  shown  in  figure  (2),  the  substrate  holder, 
dump  tube  and  endcap  were  fuzed  into  one  assembly  for  maximum 
leakage  integrity.  Not  shown  is  the  new  flat  0-ring  joints  which 
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will  replace  Che  Caper  joinCs  used  so  far. 

Up  uncil  now  rocaneters  were  used  for  main  and  bypass 
hydrogen  gas  metering  in  an  inicial  attempt  Co  keep  cosCs  low. 
These  have  been  replaced  by  elecCronic  mass  flow  eonCrollers 
which  will  give  superior  accuracy. 

The  original  plan  Co  use  Fe  or  Cr  solid-phase  doping  for 
Che  high  resiscivicy  buffer  layer  has  been  scrapped  in  favor  of 
more  promising  aeChods.  Coabinacions  of  Che  mole  fraction  effect 
and  addition  of  oxygen  during  growth  have  been  observed  to  sup¬ 
press  incorporation  of  backround  impurities  ^>7.  This  oxygen  gut¬ 
tering  occurs  without  the  aeaory  effect  caused  by  contamination 
of  Che  tube  when  heavy  metals  are  introduced  as  hi-resistivity 
dopants.  There  is  also  evidence  of  deep  oxygen  donors  further 
serving  to  reduce  the  free  carrier  concentration.  We  will  study 
other  aethods  of  introducing  oxygen  into  the  system,  e.g.,  using 
CO2  gas,  without  adversely  affecting  the  mobility  of  the  epi¬ 
taxial  active  layers. 


-  9  - 


Sumut 

We  conclude  from  the  studies  completed  thus  far  that  the 
more  conventional  designs  for  AsCl3t  which  our  redesigned  system 
more  closely  resembles,  are  best  used  for  the  growth  of  high 
purity  GaAs  buffer  and  active  layers  in  conjunction  with  "clean" 
gettering  techniques  such  as  oxygen  injection.  Various  oxygen 
compounds  such  as  CO2  and  CO  will  allow  better  oxygen  incorpora-* 
tion  and  at  the  same  time  introduce  small  quantities  of  deep  car¬ 
bon  acceptons  to  more  adequately  compensate  shallow  Si  donor  and 
other  acceptor  levels.  The  mole  fraction  effect  shows  great 
promise  in  further  dopant  suppression  when  used  in  conjunction 
with  the  oxygen  incorporation. 

With  this  considerations  taken  into  account  we  believe 
this  modified  approach  is  a  much  more  fruitful  method  to  reliably 
produce  hi-quality  epitaxial  GaAs. 
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B  Ion  I«pUBt«tion 

Semi-insulating  gallium  arsenide  grown  by  liquid  *ncip- 
suled  Cxochralski  (LCC)  method  has  been  plagued  by  a  number  of 
problems.  For  example,  heating  LEC  substrates  at  typical  anneal¬ 
ing  and  growth  temperatures  (700-900°C)  frequently  proeeduces  a 
highly  conductive  thin  layer  near  the  surface.  This  phenomenon 
is  usually  accompanied  by  poor  reproducibility  of  carrier  con- 
certrations,  low  mobilities,  and  photoluainescence  changes.  We 
have  began  to  study  this  and  other  problems  related  to  GaAs  sub¬ 
strates  for  direct  ion  implantation  applicationa .  A  simple  "get- 
tering"  technique  that  employs  heat-treatments  at  800°  for  24 
hours  in  flowing  H2  with  an  overpressure  of  As  provided  by  XnAs 
was  used.  An  outline  of  our  gettering  experimental  procedure  is 
provided  in  table  #2  on  page  15  of  this  report. 

Routine  hall  measurements  were  taken  at  room  temparature 
and  77°K.  The  hall  samples  employed  the  standard  Van  der  Pauw 
technique  with  clover-leaf  shaped  samples.  Contacts  to  the 
specimen  are  made  with  tin  beads  alloyed  to  the  leaves,  which  in 
turn  provide  edge  contacts  to  the  central  part  of  the  sample. 


For  our  purposes,  the  carrier  concentration  and  the  mobility  are 
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Figure  3b 


par ticulary  iaporCant  values,  and  can  be  obcained  from  these 
measurements . 

The  heat  treated  or  "gettered"  implanted  and  annealed 
layers  typically  showed  a  higher  room  temperature  and  77°K 
mobility.  The  improvement  in  mobility  varies  from  wafer  to  wafer 
and  is  shown  in  the  figures  (3a)  and  (3b)  of  mobility  va  sub¬ 
strate  or  boule  number.  For  instance,  Q9  showed  at  35Z  improve¬ 
ment  while  Qll  showed  no  improvement.  The  wide  range  of  data 
reveals  the  large  differences  in  the  substrate  properties, 
residual  impurities,  and  trap  concentrations. 

Standard*  photoluminescence  spectra  were  also  obtained 
for  both  the  gettered  and  ungettered  layers.  The  PL.  was  per¬ 
formed  at  7°K  using  a  Argon  layer  at  a  power  level  leas  than 
lOOmw.  The  luminescence  radiation  was  dispersed  by  a  3/4  meter 
grating  monochromator  and  detected  with  a  LN 2  cooled  photomul¬ 
tiplier  using  phase-sensitive  detection.  The  resulting  PL  spec¬ 
tra  of  the  gettered  samples  exhibited  a  reduction  in  the  Mn  and 
Cu  peak  intensities.  This  result  is  illustrated  in  the  PL  spec¬ 
trum  of  Q8  and  Q4  (figure  4). 

SIMS  analysis  experiments  were  performed  on  ungettered 
and  gettered  samples  at  Cornell  University.  The  SIMS  analysis 


c  (e-A°)  AT  7°K 


Wavelength  (A°) 


EL2 


revealed  Che  presence  of  several  residual  impurities.  These  im¬ 
purities  :  Carbon,  chromium,  boron,  managanse,  copper  and  magne¬ 
sium  were  clearly  present  in  the  10^  -  10^5  range.  In  addition, 
we  observed  outdiffusion  for  the  substrate  of  copper,  aanagenese, 
and  magnesium.  These  results  are  shown  in  figure  (5). 

Electrically  DLTS  results  figure  (6)  on  schottky  bar¬ 
riers  diodes  formed  on  the  gettered  and  ungettered  layers  reveals 
an  increase  in  the  EL2  level  and  a  reduction  in  the  broad  shallow 
peak.  We  believe  that  this  broad  peak  is  related  to  the  presence 
of  several  residual  impurities. 

In  summary,  we  have  concluded  from  the  studies  completed 
that  LEC  substrates  can  be  improved  by  heat  treating.  There 
are, however,  a  couple  of  precautions:  1)  It  is  necessary  to 
remove  all  of  the  gettered  layer  and,  2)  provide  an  suffic*nt 
overpressure  of  As.  This  work  was  presented  at  the  Electronic 
Materials  Conference  held  at  Santa  Barbara. 


GETTER  IMG  EXPERIMENTAL  PROCEDURE  TABLE  #2 

POLISH  SUBSTRATES 

DIVIDE  POLISHED  SUBSTRATES  INTO  HALVES 

BAKE  ONE  HALF  OF  THE  SUBSTRATES  AT 
800°  FOR  24  HOURS 

POLISH  BOTH  SIDES  OF  THE  GETTiERED  SUBSTRATES 
(REMOVE  35  mm) 

STANDARD  CLEAN 

IMPLANT 

Si29  6.09  X  lOl*  40  XV 

9.03  X  1011  110  KV 

4.5  X  10 12  200  KV 

ANNEAL  AT  800°C  FOR  30  MINUTES  IN  A  CAPLESS  ANNEALING  SYSTEM  WITH 
FLOWING  H 2 

COMPARE  ELECTRICAL  PROPERTIES 


!■ 

► 
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We  have  a  acquired  and  inatalled  a  MBE  system.  We  an¬ 
ticipate  that  this  system  will  be  used  to  grow  the  critical 
layers  for  our  FET  structure.  With  MBE  capabilities  it  will  be 
possible  to  investigate  normal  GaAs  FET's  as  well  as  the  high 
electron  mobility  structures  (HEMPT).  We  anticipate  that  we  will 
be  growing  both  single  as  well  as  multiple  interface  HEMPT 
devices.  Our  machine  has  the  capability  of  growing  on  2  inch 
substrates  and  we  can  load  up  to  six  wafers  in  a  single  pump 
down.  We  are  presently  baking  out  the  furnaces  in  preparation 
for  our  first  growth. 
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Figure  9 

Calibration  mark  is  one  micron 


Figure  10 

Calibration  mark  is  one  micron 


One  of  our  principal  problems  during  this  work  has  been 
the  development  of  our  lithographic  techniques.  This  section 
will  summarize  our  efforts. 

At  the  beginning  of  the  contract  a  FET  mask  set  and  a  test 
pattern  were  fabricated  at  the  National  Submicron  Center  located 
at  Cornell  University.  The  best  results  were  obtained  using  a 
wet  chemical  etch.  The  smallest  resolved  features  of  the  masks 
have  linewidths  of  better  than  0.25  microns.  We  then  went 
through  a  period  of  considerable  experimentation  with  the  expo¬ 
sure  and  development  parameters.  Our  final  process  consists  of  a 
two  layer  photoresist  technique  utilizing  a  co-polymer  of  PMMA 
and  PMAA.  The  PMMA  was  used  as  the  imaging  resist  and  develops 
more  slowly,  the  co-polymer  was  used  to  form  a  lift  off  lip 
similar  to  that  shown  in  figure  (7).  Our  final  process  is  sum¬ 
marized  in  Table  III.  Examples  of  our  better  gate  lift  off  are 
shown  in  figures  (8), (9)  and  (10). 

The  key  element  that  we  found  necessary  to  reprod uc eab 1 y 
form  .5  micron  lines  was  to  obtain  ’’good"  coformal  contact 
between  the  mask  and  the  GaAs  substrate.  In  order  to  ascertain 


whether  we  are  had  contact  a  quartz  plate  was  used.  The  substrate 
was  observed  through  a  quartz  plate  the  same  size  as  the  mask. 
Figure  (1)  is  representative  of  the  type  of  contact  we  were  able 
to  obtain.  On  the  areas  where  no  fringes  formed  we  observed  a 
darken  region  which  were  the  areas  where  we  obtained  good  con¬ 
tact.  We  found  excellent  correlation  between  these  contact  areas 
and  production  of  high  quality  lines.  As  indicated  in  Table  III 
one  of  the  process  steps  is  dissolving  the  resist  by  using  a  O2 

plasma.  For  this  process  we  have  calibrated  the  etch  rate  in  a 

\ 

barrel  plasma  reactor  and  this  is  illustrated  in  figure  (12). 


Photoresist  Etching 
Pressure  Vz  to rr 
Power  100  watt 
□  Copolymer 

a  PMMA 

o  AZ- 1350 


TABLE  III 


Proccei  Summary  for  two  layer  photoresist 


Cleaning  and  Degreasing 


Bakeout  in  Air  at  200°C  for  30  min. 


PMMA  -  PMMA  Copolymer 


is  spun  on  at  50p0A°  thickness 


Bakeout  at  160°C  for  45  min 


4Z  PMMA  is  spun  at  3000A°  thickness 


Bakeout  at  160°C  for  60  min. 


Wafer  is  placed  under  quartz  plate  to  check'  confor 


mal  contact 


Contact  adjustment 


Mask  position  checked  aad  wafer  exposed  for  6  minutes  at 
a  lamp  intensity  of  19.8aw/cm2 
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Sample  is  developed  in  solution  of  2  parts  Propanol  1 
part  Toluene 

Oxygen  Plasma  Clean~up  1  minute  at  100  watts  a  1/2  torr 
pressure 

Sample  is  now  ready  for  metallization 


Device  Structure  and  Results 


In  this  section,  we  will  summarizies  our  device  results. 
Since  as  mentioned  in  a  previous  section  our  vapor  systems  were 
not  fully  operational  all  of  devices  results  are  on  ion  implanted 
material.  This  material  was  baked  out  in  order  to  attempt  to 
getter  residual  impurities. 

Our  first  submicron  structure  attempted  to  illustrated  our 
lithography  technology  using  a  self  aligned  Ti/*.-  gate.  In  this 
structure  the  gate  was  deposited  first  and  the  source  and  drain 
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were  aligned  around  the  gate.  This  structure  was  alloyed  with 
the  gate  in  place.  We  found  it  necessary  to  deposited  a  minimum 
of  2000”A  of  Ti,  to  provide  adequate  protection  during  our 
source/drain  alloying  process.  In  the  future,  we  intend  to  ex¬ 
periment  with  Pt  as  a  barrier  metal.  Using  this  fabrication 
sequence  we  have  fabricated  working  devices.  Data  and  pictures 
for  representation  devices  are  shown  in  Table  IV  and  figures 
(  13)  ,(  14)  ,and  (15)  respectfully.  We  observed  some  problems 
with  edge  effects  on  the  ohmic  contact  metalization  layers.  We 
beliewe  these  problems  can  be  corrected  by  the  addition  of  nickel 


to  the  metalization  system 


Figure  14 


i/division  2  ma/division 


Figure  15 

I-V  characteristics  of  .5  micron  FET  (top) 
I— V  characteristics  of  gate  of  FET  (bottom) 


After  demonstr ating  working  devices  we  proceeded  to 
develop  a  standard  process  for  our  short  gate  FET  experiments. 

In  this  structure  again  the  gate  is  deposited  first  then  a  "T" 
structure  is  foraed  by  selectively  etching  the  Ti  in  a  CF4  plas¬ 
ma.  This  is  followed  by  a  blanked  evaporation  of  Au-Ge. 

The  source  drain  regions  are  then  defined  aaking  this 
structure  a  self  aligned  gate  as  well  as  self  aligned  source 
drain  structure.  We  antic ipate ■ that  this  structure  should  lower 
the  parasitic  resistance  in  the  source  region  enough  to  enable 
isolation  of  the  effects  of  gate  length  on  device  performance. 
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Theorectical  Studies 


We  have  investigated  several  areas  related  to  PET  device 
performance.  The  details  of  these  investigations  appear  in  the 
following  papers. 
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Abstract 


Whereas  improvement  of  the  interface  between  the  active  layer  and 
the  buffer  layer  has  been  demonstrated,  the  leakage  effects  can  be 
important  if  the  buffer  layer  resistivity  is  not  sufficiently  high 
and/or  the  buffer  layer  thickness  is  not  sufficiently  small.  We  find 
two  buffer  leakage  currents  exist  from  the  channel  under  the  gate 
to  the  source  and  from  drain  to  the  channel  in  addition  to  the  buffer 
leakage  resistance  between  drain  and  source.  It  is  shown  that  for  1  p 
gate-length  n-GaAs  MESFET  if  the  buffer  layer  resistivity  is  12«a-*Bm 
and  the  buffer  layer  thickness  h  is  2  p,  the  performance  of  the  device 
degrades  drastically.  We  suggest  that  h  should  be  below  2  ^a. 


I .  Introduction 


The  main  advantage  of  using  a  buffer  layer  in  n-GaAs  MESFET  is 
that  the  interface  between  the  active  layer  and  substrate  can  be 
improved.  However,  since  the  buffer  layer  placed  between  the  active 
layer  and  the  substrate  has  'thickness  of  several  microns  and  resistivity 
much  higher  than  that  of  the  active  layer  and  lower  than  that  of  the 
substrate,  leakage  currents  may  flow  from  and  to  the  active  layer 
through  the  buffer  layer  due  to  its  lower  resistivity.  The  leakage 
effects  on  dfcvice  parameters  and  in  particular  on  device  performance 
will  be  examined  in  this  work. 

In  previous  works  dealing  with  leakage  in  a  semi-infinite  substrate 
(liechti,  1976  and  Reiser,  1973)  the  model  of  the  leakage  current 
was  considered  simply  as  a  leakage  resistance  between  the  source  and 
drain.  Whereas  this  simple  model  is  useful  for  explaining  the  effect 
of  the  leakage  resistance  on  the  1-7  characteristics,  a  more  detailed 
model  is  desirable.  In  this  work  we  have  considered  the  leakage 
problem  of  a  finitely  thick  buffer  layer  and  made  detailed  calculations 
which  enable  us  to  present  a  better  model.  In  addition  to  the  drain- 
to-SQucna  .leakage  resistance,  the  channel-to-source  and  drain-to- 
channel  leakage  currents  are  included  in  our  model.  Furthermore  we 
have  considered  in  our  model,  the  varying  channel  potential  which 
is  of  importance  for  short  gate  FET. 

In  order  to  simplify  the  calculation,  we  have  neglected  the 
possible  effects  of  interface  trapping  and  the  space  charge.  Also 
we  assume  that  the  diffusion  current  from  the  active  layer  to  the 
buff er  layer  is  negligible  and  that  the  transverse  field  in  the  channel 
is  much  smaller  than  the  longitudinal  field  so  that  the  voltage 
along  y-direction(transverse  direction) can  be  considered  to  be  the 


same  across  the  channel  -  This  implies  that  the  voltage  on  the  channel- 
gate  interface  is  the  same  as  that  on  the  channel-buffer  layer 
interface  under  the  gate.  Thus  a  virtual  gate  electrode(see  Pig.  1) 
can  be  assumed  to  exist  along  the  interface  between  the  active  and 
buffer  layer.  In  addition  it  is  assumed  that  the  gate-drain  and  gate- 
source  spacing  are  negligibly  small.  Based  upon  these  assumptions 
it  is  therefore  possible  to  deal  with  merely  a  static  field  problem 
in  the  buffer  layer.  Along  the  interface  between  the  buffer  layer  and 
S. I. substrate,  no  leakage  is  assumed  to  exist,  i.e.  the  substrate  is 
assumed  to  be  a  perfect  insulator;  leakage  occurs  only  in  the  buffer 
layer  and  its  interface  with  the  active  layer. 

By  solving  the  field  problem  with  the  help  of  conformal  mapping 

we  can  calculate  the  buffer  leakage  resistance  rds^  and  the  two 

leakage  current  sources  1^,  and  1^,  on  using  superposition  principle 

for  current  on  the  three  electrodes(Sec.  II).  With  this  model  in 

mind  we  derive  the  equivalent  circuit  including  the  leakage  effect 

0 

from  which  the  h-parameters  for  the  device  with  and  without  buffer  layer 
can  be  determined.  The  performance  parameters  such  as  gains  and  PMty 
can  be  obtained  and  compared  for  MESPET's  with  and  without  buffer 
layer (Sec.  III).  Numerical  results  and  conclusions  are  presented  in 
sections  IV  and  V*  respectively. 


II.  Analysis  of  Leakage  Elements 


We  will  derive  the  leakage  elements ,  l.e.,  leakage  resistance  and 
the  leakage  current  sources,  hy  solving  two-dimensional  Laplace 
equation  for  potential  distribution  in  the  buffer  layer  with  appropriate 
boundary  conditions  set  on  the  interface  between  the  active  and 
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buffer  layers.  In  Fig.1,  S  and  L  are  the  virtual  source  and  drain 
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electrodes  and  G  is  the  virtual  gate  electrode  on  which  the  voltage 
varies  along  x-direction  due  to  voltage  variation  in  the  channel. 

Since  the  field  solution  assiciated  with  Fig.1 i with  finite  buffer 
layer  thickness  h  is  difficult,  we  will  use  conformal  mapping  three 
times  to  transform  this  original  structure  (A)  to  the  parallel  plate 
structure  (c)  as  shown  in  Fig.2-‘c  from  which  the  field  and  potential  . 
can  easily  be  obtained. 

As  shown  in  Fig.2-a  (structure  A)  S  ,  G  and  D  are  the  virtual 
source,  gate  and  drain,  respectively  where  source  and  drain  are 
extended  into  infinity  along  ±x  and -the  buffer  layer  has  a  finite 
thickness  h.  Structure  A  can  be  conformally  mapped  into  structure  B* 
(Overmeyer,  1970)  with  the  transformed  dimensions  given  by 


*2*  -^tanh(-^---) 


c  -  £— 

c2  TT 


(D 

(2) 


where  a^  is  one  half  of  the  given  gate  length.  A  further  transformation 
gives  structure  B  with  new  dimensions  as  given  by 


or 
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x 

o  o 


7t  ai 

tanh(  2%  J 


(3) 


(4) 
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Prom  structure  B  the  final  parallel-plate  structure  (c)  is  obtained 
again  by  using  conformal  mapping  (Okoro,  1980).  Prom  structure  c  we  have 

L®2K(p)  (gate  length)  (5) 

'  T»K* (p)  (buffer  layer  thickness)  (6) 

where  K(p)  and  K?(p)  are  the  elliptical  integral  of  the  first  kind 
and  p  is  given  in  (4). 

When  bias  voltages  are  applied  to  the  gate  and  drain,  the  channel 
potential  will  vary  along  the  longitudinal  direction(from  S  to  D  ). 

As  a  result  of  conformal  mapping  the  field  variation  is  preserved, in 
the  various  structures.  The  virtual  gate  potential  distribution  in 
the  original  structure (structure  A)  can  be  approximated  by  two  linear 
regions  based  on  the  model  of  Pucel,  1975(see  also  Wang,  1979)  as 
shown  in  Pig.  3  where  7q(0)*Q,  VQ(<x)=7p  and  7Q(2a1  )*3TDS  where  7^  and  - 
VDS  are  the  pinch-off  and  drain  to  source  potentials,  respectively,* 
is  the  position  of  pinch-off  point  and  2a1  is  the  total  given  gate 
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length.  The  virtual  gate  potential  distribution  on  structure  C  can 
be  transformed,  from  Pig.  3  into  Pig.  4  toy  requiring  that  the  source, 
pinch-off  and  drain  voltages  remain  invariant, i.e.  Vq(0)»0,  7q(*)» 7^ 
and  70(2K(p))»7^s  Where  I«2K0p)  is  the  length  of  the  parallel-plate 
structure  (C).  In  region  I  of  Pig.  3  we  have 

7Q(x)»m1x  (7) 

where  m^P^/oc  and  in  region  I  of  Pig.  4  we  have 

* 

7o4±)»m^x  (8) 

where  .  But 

oc  ok 


(9) 


(10) 


from  which  m1  can  he  obtained,  and  is  given  by 


Similarly  in  region  II  of  Pig. 3  we  obtains 

V^xJ-m^x-cO+Vp  <*4x^2ai  (11) 

where  m2“(VDS-Vp)/(2a1- ,  and  in  region  II  of  Pig. 4  we  obtain 

Y§(±)«m^(x-©c)+Vp  2K(p)  (12) 

where 

,  2a.,-* 

m2*m2  SKip'J-o*' 


(13) 


The  parameters  for  the  two  linear-region  approximation  of  gate  . 
potential  distribution  are  summarized  in  table  1 . 

Here  we  will  make  use  of  superposition  principle*  by  which  we  will 
block  one  of  the  three  electrodes, i^e. ,  source,  gate,  or  drain,  at  a 
time  and  calculate  the  leakage  element  so  that  the  three  leakage 
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elements, i.e. ,  rdal)  ,  I^g  and  Ij^,,  can  be  determined  from  the  final 
parallel-plate  structure. 

A.  leakage  resistance  rdafe 

Leakage  resistance  is  defined  as  the  resistance  between  source 
and  drain  when  the  gate  electrode  is  blocked,!. e. ,  there  ia  no 
current  flow  on  the  gate  or  (3V/ay) jy,^*0  as  shown  Pig. 5.  Therefore, 
the  structure  reduces  to  a  parallel-plate  structure  in  which  the 
leakage  resistance  between  drain  and  source  is  given  by 
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where  f  is  the  resistivity  of  the  buffer  layer ,  2K(p)  is  the  distance 
between  drain  and  source  and  A  is  the  cross-sectional  area  which  is 


given  by 


A-K  (p)z 


(15) 


where  z  is  the  width  of  the  structure  (in  z  direction).  Thus  the 
leakage  resistance  in  (14)  becomes 


*dsb  z 


t  2K(P) 

*  tfp) 


(16) 


B.  leakage  currents 


Leakage  between  the  gate  and  source,  and  between  the  drain  and 

i  t 

gate  are  modelled  as  two  leakage  sources  1^  and  instead  of  two 

resistors,  because  on  the  gater  the  voltage  varies  along  x-direction 
and  therefore  uniform  leakage  resistances  can  not  be  defined.  On  using 

i 

the  superposition  principle  IGS  is  obtained  by  calculating  the  voltage 
distribution  in  the  structure  C  with  the  drain  electrode  being  blocked, 
and  is  obtained  with  source  blocked. 

(i).  Source-to-gate  leakage  current  1*^ 

The  structure  C  when  the  drain  is  blocked  is  shown  in  Big.  6.  Note 
that  there  is  no  leakage  between  the  buffer  and  the  S-I  substrate. 

The  general  eigen  solution  of  Laplace  equation  as  applied  to  Fig. 6 
is  given  by 

Vn( x , y ) « £aqc os hC^y^+B^ inh (*ny  )j  ^Cacos (<*nx)+Dnsin(<*nx)j  (17) 

Using  the  boundary  conditions  shown  in  Fig.6,i.e., 

v*,y)L0'0  (is) 


(18) 


where  7Q(x)  is  the  channel  and  therefore  the  buffer  layer  potential. 
The  potential  distribution  in  the  buffer  layer  can  be  obtained  and  is 
given  by 


v(x,y) 


sin  cosh  2*Z 


where 


~ZWJ  £odd  ccsh[Sf-$$] ~  ««P> 

?K(p) 

n  a  j  7  Main  h'X.  X  _  dX  (23) 

n  X  0 


IK(p) 

,(x)sin  -n7r-X-  43 
'b  J  4K(p)  *  . 


or  from  Fig.4- 


£K(p) 


4K(p) 


(24) 


Substituting  the  appropriate  values  of  Y«(x)  as  shown  in  (8)-and  (12) 
into  (24)  •  we  obtain 

Z,n2sin  T£+  [^i]2(nl''l,,2)sin  TOlT  (25) 

t 

The  total  leakage  current  I^g  is  given  by 


T*  z 
XGS  7 


Mv) 


iZlj 


dy 


(26) 


x»0 


where  z  is  the  width  of  the  structure  and  ,/  is  the  resistivity  of  the 


buffer  layer.  On  using  (22),  (25)  and  (26)  the  total  leakage  current 


from  the  gate  to  source  can  be  obtained  and  is  given  by 


^K(p)  n«odd 


untanh[sjl 


nr  «p) 
^  K(p). 


(27) 


(ii).  Drain-t o-gate  leakage  current  I. 


The  structure  C  when  source  is  blocked  is  shown  in  Fig. 7.  Fig.7-a 
shown  the  drain  is  biased  at  tds-  This  is  equivalent  to  the  situation 
as  shown  in  Fig.  7-b  from  which  a  solution  similar  to  that  in  the 
previous  subsection  canrbe  obtained.  Thus  using  boundary  conditions  as 


shown  in  Fig. 7-b,  we  have 


K(p).V(x,y)» 


n-odd 


cosh 


[ill 


_ mf  x  _ nlfy  (28) 

008  WpT  008 


where  un  is  given  by 


VSi+,4ii'J  [Vx>-7ds]  008 


(29) 


On  using  the  appropriate  values  of  7Q(x)  given  in  (8)  and  (72)  in 
(29) ,  we  obtain 


vt4^]2^-®^008  TOify  *®i] 


(30) 


And  the  total  drain  to  gate  leakage  current  1^  is  given  by 

.  J*(P) 

h»m  77~J0'  ^f*7*  d7 

0  x-2Z(p) 


(31) 


hoT  'TTpT 


n»odd 


u^sin  S^tanh^S^L  ^ 


(32) 


Note  that  the  total  leakage  currents  given  by  (27)  and.  (32)  are 
inversely  proportional  to  the  buffer  layer  resistivity  which  implies 
that  higher  buffer  layer  resistivity  will  give  rise  to  a  decrease 
in  leakage  currents  resulting  in'  better  device  performance.  However, 
as  will  be  shown  later,  the  buffer  layer  thickness  is  critical  to 
leakage  resistance  and  currents. 

It  should  also  be  noted  that  the  gate-to-source  and  the  drain- 
to-gate  leakage  currents  are  the  same  as  the  channel- to-source  and 
drain-to-channel  leakage  currents,  respectively. 


III.  Performance  Parameters 


The  equivalent  circuit  of  MESPET  including  the  leakage  elements 
obtained  in  Sec. II  is  shown  in  Pig. 8.  The  performance  parameters  of 
the  device  can  be  obtained  using  h-parameters(Ohkawa,  1975).  Using 
simple  network  theory,  the  h-parameters  of  Pig. 8  can  be  obtained^-as  follov 


*7 


Y„«0 


g3fcj _ 

>=g3^"=dg-“Cg3ri°dg 


(33) 


i2  ~Tl 


5df 


I^O 


R3 


**21' 


V° 


gB-3"gdgO*>»eg3gl> 

*'egs+J"'Wi+*'eg«riJ 


(3*) 


(35) 


h  i 
Pi*? 


(°dg/C6S:‘B)Ug»>':!"C<ig) 


+(: 


'dbs 


,+^cdg)  (36) 


Where 


A« 


SgB~S<Sl 


J«c 


f 


gs 


*dbs  ?dsb+3wCdsrdsbRds+ilds 


®dsrldsb 


(37) 

(38) 

(39) 


Por  the  intrinsic  device  parameters  see  Table  2.  The  drain  to  gate 
leakage  trans conductance  gdg  is  given  by 


11 


gdg" 


lDG 


'DS 


(40) 


VGS  ■constant 


When  no  leakage  exists  in  the  buffer  layer,  1^,  and  I^g  are  equal  to 
zero  and  *> ,  the  resultant  h-parameters  approaches  to  the  well- 

known  intrinsic  h-parameters (Ohkawa,  1975). 

The  performance  parameters  are: 

A.  Maximum  stable  gain  GMS 


GMS* 


|h12| 


(41) 


B.  Unilateral  gain  U 


h12+h2ll 


4 [lie  t tLy  t  i-fte  ( hg 2 )  +Tm ( K[ 2 ) Im { A  )] 


(42) 


where  He  denotes  the  real  part  and  Im  denotes  the  imaginary  part. 
G-.  Maximum  frequency  of  oscillation  FMAY 


MAX 


2VrAis' 


<«) 


17.  Numerical  Results 


The  MESFET  numerically  analyzed  in  this  section  has  a  gate 
length  of  1  jx  and  a  width  of  500  ji  with  an  active  layer  thickness 
of  0.2ju  arid  a  buffer  layer  resistivity  of  12 -ft— cm.  The  typical 
intrinsic  parameter  values  for  the  MESFET  is  shown  in  table  2 
(Liechti,  1976).  Unless  otherwise  stated,  the  buffer  layer  thickness 
is  assumed  to  be  2  jx. 

The  comparisons  between  the  performance  parameters  of  the  MESFET 
with  and  without  buffer  layer  are  based  on  the  intrinsic  equivalent 
circuit  in  which  the  parasitic  elements  have  been  neglected(see  Fig.8) 

The  leakage  current  sources  are  strongly  dependent  on  the 
resistivity  of  the  buffer  layer  as  shown  in  (27)  and  (32).  Therefore % 
by  increasing  the  .resistivity  of  buffer  layer,  the  leakage  currents 
will  decrease  which  improves  the  device  performance.  In  Fig.9  the 
leakage  resistance  and  currents  are  given  as  functions  of  the  buffer 
layer  thickness.  In  order  to  reduce'  the  leakage,  the  buffer  layer 
thickness  must  be  less  than  2  jut.  If  otherwise  leakage  will  short-circui- 

t  » 

the  active  layer  current.  It  is  noted  that  I&g.  This  can  be 

accounted  for  by  the  existance  of  high  field  region  between  the  drain 
and  gate.  The  variation  of  the  leakage  currents  with  the  normalized 
drain  bias  voltage* ^DSN  is  silown  in  Jig.  10,  where  the  slopes  of  the 
curves  give  the  leakage  transconductances. 

It  is  shown  in  Fig. 11  that  F^A^  for  the  buffered  device  degrades 
by  a  factor  of  2  to  4  due  to  leakage  effect.  By  increasing  the  buffer 
layer  resistivity  and/or  decreasing  the  buffer  layer  thickness 
will  improve  proportionally. 


VI.  Conclusions 


The  leakage  effect  in  MESFET  has  been  investigated.  It  has  been 
shown  that  if  the  buffer  layer  resistivity  is  low  and/or  the  buffer 
layer  thickness  is  large ,  leakage  effect  will  considerably  degrade  the 
device  performance.  We  have  shown  that  for  the  1  ji  gate-length  n-GaAs 
MESEET  if  the  buffer  layer  resistivity/  is  12A-cm  and  the  buffer 
layer  thickness  h  is  2  ja,  the  gains  of  device  decrease  drastically 
from  their  intrinsic  values.  Since  the  leakage  resistance  is  proportior 
to/  and  the  leakage  currents  are  inversely  proportional  to/  thus 
leakage  effects  can  be  reduced  by  increasing/.  However,  both  the 
leakage  resistance  and  currents  depend  on  h  more  strongly  than  /  thus 
it  is  advantageous  to  decrease  h  for  improving  performance.  We 
suggest  that  h  should  be  lower  than  2  /. 


Table.  1 


Parameters  for  two-linear  region  approximation  of  gate  potential 

distribution 


(known; 

Original  structure 

- (calculated) — 

Pinal  structure 
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m1 

»  Tp 

°1  " 

“2 

»  Do  ^0 

nu-  — - 7— 

2K(p)-<* 

Table.  2 

Typical  intrinsic  parameter  values  for  the  MESFET  with  the  gate  length 
of  1  jx9  width  of  500  ji,  active  layer  thickness  of  0^2  with  the  doping 

x  .  ..  ,  density  of  lO1*^  cm“^ 
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1.  Assumed  MESFET  structure  where  S  ,  G  and  D  are  the  virtual 
source,  gate  and  drain,  respectively* 

2.  Conformally  mapped  'structures.  Final  structure  is  the  parallel- 
plate  (c). 

3.  Potential  distribution  of  structure  A* 

A.  Potential  distribution  of  structure  C. 

3*  Leakage  resistance  between  source  and  drain  when  gate  is  blocked 
(non-conducting). 

6.  Leakage  current  from  gate  to  source  when  drain  is  blocked. 

7.  Leakage  current  from  drain  to  gate  when  source  is  blocked. 

8.  The  ac  equivalent  circuit  of  a  MESFET  with  leakage  in  the  buffer 
layer,  where  IDG-gdg.72  1GSmgg3m72* 

9*  Plot  of  leakage  elements  vs.  h  for  1  jx  gate  length  MISFET  when. 
V^g— 1.44  volts. 

1  (X< . Total  leakage  currents  vs.  for  1  jx  gate  length  MESFET  with 

buffer  layer  thickness  of  2  jx  and  buffer  layer  resistivity  12-a--cm 

1 1 .  Fmay  vs.  V^gjj  for  MESFET  with  t  jx  gate  length  and  the  buffer  layer 
with  the  resistivity  of  12A-cm  and  the  thickness  of  2  jx. 

12.  Unilateral  gain  vs.  frequency  for  1  /i  gate  length  MESFET  with  the 
buffer  layer  thickness  of  2 jx  and  resistivity  of  12J*--cm. 

13.  Maximum  stable  gain  vs.  frequency  for  1  jx  gate  length  MESFET  with 
the  buffer  layer  thickness  of  Z  ix  and  resistivity  of  12-a  -cm. 
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ABSTRACT 


An  analytical  model  of  the  GaAs  MESFET  with  arbitrary 
non-uniform  doping  is  presented.  Numerical  results  for  linear 
lateral  doping  profile  are  given  as  a  special  case.  Theoretical 
considerations  predict  that  better  device  linearity  and  improved 
can  be  obtained  by  using  linear  lateral  doping  when  doping 
density  increases  from  source  to  drain. 


I.  INTRODUCTION 


Improved  linearity  and  noise  figure  in  FET's  with  vertical 
non-uniform  doping  has  been  demonstrated,  by  Williams  and  Shaw(1973), 
Roberts,  Lynch,  Tan  and  Gladstone (1 973) ,  Pucel(1978)  and  Santis (1979). 
In  this  paper  we  consider  the  non-uniform  doping  (along  the  electron 
motion)  in  the  GaAs  layer  of  the  PST.  Two  advantages  in  this  type 
of  PET  are  anticipated,  First  if  an  appropriate  profile  is  U3ed, 
the  depletion-channel  interface  can  become  more  uniform  along  the 
lateral  direction  hence  better  linearity.  Secondly,  The  cut-off 
frequency  P^  increase  can  be  realized  for  certain  lateral  doping 
profile. 

A  general  mathematical  analysis  is  presented  in  which  the 
depletion  region  potential  is  solved  from  Poisson's  equation  as 
a  boundary  value  problem.  In  the  channel  we  adopt  the  saturation 
velocity  model  introduced  by  Williams  and  Shaw(1978)  which  is 
valid  for  short  gate  (Sze,  1981).  The  analysis  takes  into  account 
two-dimensional  general  non-  uniform  doping  profile.  In  obtaining 
numerical  results  however,  a  simple  linear  doping  profile  along 
the  lateral  direction  is  assumed.  Device  parameters  such  as  gm, 

Cgg  and  P,j  are  given  in  terms  of  the  gate  bias  voltage  for  the 
PET  with  gate  length  of  0.6  jx* 

In  Sec.  II  the  two  dimensional  potential  in  the  depletion 
region  is  analyzed.  In  Sec.  Ill  the  depletion  height  h/a,  reduced 
potential  and  channel  current  for  linear  lateral  doping  case  are 
presented.  Device  parameters  are  obtained  in  Sec.  IV.  Conclusions 
are  given  in  Sec.  V. 


z 


II.  TWO  DIMENSIONAL  ANALYSIS  OP  THE 
DEPLETION  REGION  POTENTIAL 

A.  General  solution 

The  two  dimensional  Poisson's  equation  is  given  by 

V2V(x,y).  a2y<ftr?  +  -f(x.y)  (1) 

Si  3y 

where  f  (x,y)*qN(x,y)/«  ,  N<x,y)  is  the  x-and-y  dependent  doping 

density  and  W(x#y)  is  the  potential  in  the  depletion  region  as 
shown  in  Pig.  1 .  This  inhomogeneous  second  order  partial  differential 
equation  can  be  solved  by  using  Green's  function  method  where 
equation  ( 1 )  will  be  written  as 

«-a(x-r)a(y-’>)  (2) 

s 

where  f  and  rj  are  the  dummy  variables.  The  boundary  conditions 
associated  with  (1)  are  given  as  follow 


Schottky  gate: 

W(x,y) 

-0 

y»0 

(3) 

Depletion-channel 

<♦} 

Interface : 

ay 

-0 

y«h 

Ohmic  contact  source: 

ax 

-0 

(5) 

x»0 

Ohmic  contact  drain: 

S¥(x,y) 

dx 

H 

■ 

*** 

u 

o 

(6) 

where  h  is  the  variable  depletion  layer  height  and  /  is  the  gate 
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length  (see  Jig.  1).  It  is  the  unique  property  of  the  operator 
(v2)  which  is  the  sum  of  the  two  commutative  operators  which  enable 
ua  to  find  an  explicit  solution  to  (1).  Equation  (1)  can  be  written  as 

72W(x,y)»K1W(x,y)+K2W(x,y)  (7) 

where  K1*  d2/dx2  and  Kg*  d2/dy2.  Friedman  (1956)  shows  that  if  K1 
and  Kg  commute,  Kg  (or  K^)  can  be  treated  as  a  constant  so  that  the 
partial  differential  equation  (1)  becomes  an  ordinary  differential 
equation  given  by 

jO.fei.yl  +  m2W(x,y)«- f(x,y)  (8) 

.  dx2 

where  m»  ^K^.  Using  the  boundary  conditions  given  in  (5)  and  (6), 
the  solution  of  (8)  becomes 

w(xfy)«  iy*-f(r»y)aih[m(x-r)]dr  (9) 

where  f(f,y)  is  the  spectral  representation  of  Kg.  The  eigenvalues 

of  will  be  purely  imaginary  i,e.f  m*(3n»)/2h,  and  consequently 

the  sin  term  in  (9)  will  become  sinh.  The  eigenfunction  of  K1  can 
be  determined  by  considering  the  following  homogeneous  ordinary 
differential  equation 

-d.*{xrY)  +x2w(x,y)»0  (10) 

dy 

where  .  *VK1.  Using  boundary  conditions  given  in  (3)  and  (4),  the 
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eigenfunctions  of  W  can  be  obtained  and  they  would  be  of  the  form 
sinxy  where  x«(n*y)/2h  (n»1 ,3,5, . . ) .  Because  of  the  fact  that  W 
is  orthogonal  and  it  can  be  normalized  to  unity,  let 

«r.y)-Z  [/5n8lnS2|-l  (11) 

n«oddL  J 


where 


(12) 


where  2/h  is  the  normalization  factor.  Substituting  (12)  into  (11), 
the  spectral  representation  of  i ^  can  be  obtained  and  is  given  by 

1 S  8 

nsodd 


Note  that  f(f,>?)  includes  the  general  doping  density  in  dummy 
variables  f,i»  ,  i.e.,  f  (x,y)»qN(x,y)/«  .  The  depletion  region 
potential  W(x,y)  of  the  PET  with  general  two  dimensional  doping 
density  can  be  obtained  by  substituting  (13)  into  (9)  and  it  is 
given  by 


W(x,y): 


s 


(1/n)sin(||Z)  |y*dr  N(r,n) 
in(§£2)sinh[§£  (x-r)]j 


ia  V 

r«  /  j 
Ha  Odd 


(H) 
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B.  Uniform  doping  profile 

Whereas  equation  (14)  satisfies  (1)  and  its  boundary  conditions 
as  can  be  verified,  it  reduces  to  the  solution  for  uniform  doping 
(Pucel,  1975)  by  setting  N(f,?)-N0  in  (14).  Thus  at  the  interface 
(y**h)  equation  (14)  becomes 

W(x,h)=-  ~  V]  (l/n)sin(  -*p)  f  f iv  sin(  ^^)  a  inhC^C*”  HI  (15) 

n«odd  J  J0 


After  integration  and  some  algebra,  equation  (15)  becomes 


16qh^Nn  r — r  , 

W(x,h)= — — 2^  (1/a5)sin(-^-) 


n«odd 


The  identity 

^*  (  1  /n^ )  sin(  -^)  ■  yg 
n>odd 


(16) 


(17) 


can  be  applied  to  (16)  to  give  the  depletion  region  potential  in 
a  uniformly  doped  PET,  i.e.. 


w(x,h).w00(  |  )2 

which  is  identical  to  that  obtained  by  Pucel  (1975). 
pinch-off  potential  of  a  uniformly  doped  PET. 


(18) 

*00  18  the 

(19) 


Note  that  there  is  no  well-defined  pinch-off  potential  for  non- 
unif ormly. doped  PET. 


III.  LINEAR  LATERAL  LOPING  PROPILE 


For  the  special  case  of  lateral  doping,  N(f,»?)»N(f)  in 
equation  (14),  where  N(0  is  the  lateral  doping  profile.  Therefore, 
the  depletion  region  potential  for  lateral  doping  becomes 

x  -h 


W(x,y)»-  ^  (1/n)sin(  |j£)  f d?f dn  N(f)sin(  ||^)sinh[|5(x-r)]  (20) 

n  « odd  ® 


A.  Lepletion  height  and  reduced  potential  distribution 
For  linear  doping  let 


N(O»N0(1  +  «$) 


(21) 


where  *  is  the  rate  of  change  of  doping  density  which  increases 
with  f  (or  x)  if  «>o  and  decreases  if  «<0.  Nq  is  assumed  to  be 
1017  cm”5.  Substituting  (21)  into  (20), 

Vf(x,y)~  (1/n)sin(  f^)  f At  fit  (1  +  «?)sin(igl)3lnh[§§(x-S)] 

nsodd  (22) 


After  integration  with  respect  to  ,  (22)  becomes 

8qhN, 

1  / 
m 

nsodd 


W(x,y)»- ^2  (1/n2)sln(  ffZKA^A-,)  (25) 

^  ^  _ I  ■ 


where 


A1*  sinh[g£(x-/X! 

Ag-  8inht^x-l )] 


(24) 

(25) 
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Thus 

A  gh 

Ai  mr 

(28) 

and 

•  2h*x 

a2  5 W 

(29) 

Therefore  equation  (23)  becomes 
8qhN0 

ns  odd 
or 


v(*.y)~  ^2  V  (1/n2)ain(ffE)<- 
7T  £  JrrL, 


2h  _  2h«x 
n?f  “  mr 


) 


(30) 


16qhZN,v 

WU.y) - J-2  >,  (1/n3)sm(  §gZ)(1-H«)  (31) 

K  £  nsodd 

At  depletion-channel  interface,  i.e.,  y»h,  the  depletion  region 
potential  W(x,y)  becomes 

W(xfh)«  — 3~H-(1  +«x)  2^  (l/n3)sin(-f2T)  (32) 


"sodd 


Using  the  identity  described  in  (17) »  equation  (32)  reduces  to 

qh2N, 


W(x,h)«  -gg-y.(i+ctx) 


(33) 


The  reduced  potential  u  which  is  the  same  as  w  of  Pucel  (1975)  can 
be  obtained  as  follows 


u2(*,h).  'OgSl  .r^l2(l«a) 


oo  L  a  J 


<34) 


where  WqQ  is  given  in  (19).  Because  Splays  an  important  role  in  our 
analysis,  its  range  and  limitation  need  to  be  considered.  When 


otis  too  small,  the  doping  becomes  almost  uniform,  however,  o<  can 
not  be  too  large  for  otherwise  the  semiconductor  will  become 
degenerate.  We  shall  consider  the  two  cases, i.e.,  small  *0.1) 

and  large  c*  ( J o< J^j  =0, 9^  where  the  gate  length  is  assumed  to  be  0.6/4 

(i) .  Small  <X(j©t^|  =0.1) 

In  this  case  the  normalized  depletion  height  h/a  becomes 
almost  uniform  in  the  saturation  velocity  model  (Sze,  1981), i.e., 
u=s  where  s*u(0,h).  Therefore  (34)  becomes 

*  s  (35) 

Where  s  is  the  reduced  potential  at  the  source  and  is  given  by 
(Pucel,  1975) 

Where  Vgg  is  the  gate-source  bias  voltage  and  <P  is  the  barrier 
potential. 

(ii) .  Large  c<  (JotJ?J»0.9) 

In  this  case  It  Is  necessary  to  calculate  the  potential 
V(x)  in  the  channel.  Prom  Pig.  1  and.  equation  (34), 

v(x)»  -(Vgs+ <<>)+w(x,h)»W00[u2(x,h)-s2]  (37) 

Under  the  commonly  used  assumption  of  neutral  channel,  where  the 

carrier  and  doping  densities  are  the  same,  the  Poisson’s  equation 
becomes 


The  assumption  of  neutral  channel  is  reasonable  for  high  doping 
density.  It  is  obvious  that  the  solution  of  (38)  is  V(x)«Ax-*-3, 


Where  A  and  B  are  constant  to  be  determined.  Using  the  conditions 
that  V(0)*0  (source  is  grounded)  and  V(tf)»Vds  (drain-source  bias 
potential), .one  obtains  •  B-O,  -  and  A=Vds/i.  Therefore  the  channel 

potential  V(x)  becomes 

V(x)*(Vds/J?)x  (39) 


Prom  (37)  and  (39), the  reduced  potential  u(x,h)  can  be  obtained 
and  is  given  by 

u(x,h) 

And  from  (34)  the  normalized  depletion  height  h/a  can  be  obtained 


•/ 


ds7 


rOO 


(40) 


and  is  given  by 


a 


T+c[x 


(41) 


Where  7d3n*Vd3/W(X).  Equation  (41)  implies  that  as  <X  increases,  the 
height  of  the  depletion  region  is  no  longer  a  constant  and  varies 


with  <*  and  x. 

B. Channel  Current 

The  total  current  consists  of  the  conduction  current  and 
diffusion  current.  However,  in  the  neutral  channel  the  diffusion 
current  can  be  neglected.  Also  neglected  here  is  the  small  band- 
gap  narrowing  effect  due  to  doping  variation. 

The  conduction  current  is  given  by 


I(x)-qn(x)vsA 


(42) 


-t 


Where  q  is  the  electron  charge,  n(x)*UQ(1  +  «x)  is  the  carrier 
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L  (or  doping)  density,  vs  is  the  saturation  velocity  and 
Q  A(x)*z[a-h(x)]  is  the  cross-sectional  area(z  is  the  device  width), 
|:  Thus  from  (42)  we  have  . 


I(x)-qN0(1+  Ux)  [l-  azvs  (43-a) 

In  equation  (43-a) (which  is  equivalent  to  (6)  of  Pucel,  1975) »  the 
current  is  a  function  of  h(x)  and  n(x).  Therefore  the  total  average 
current  I  ('which  is  constant)  is  obtained  by  integrating  (43-aJ 
from  x«0  to 

(43-b) 

Thus  for  small  (|<rtJ?|=*0.1),  h(x)/a«s  (see  (35)),  and  the  total  average 
current  becomes 

I«qvflzaN0(1-s)(1+  %£-)  (44) 

4 

And  for  large  of.  (jot#  J»0.9) ,  the  total  average  current  can  be  obtained 
from  (43)  »  and  is  given  by  ^ 

I  -  SM]  dx  (45) 

It  is  convenient  that  we  preform  the  integration  in  u  instead 
of  x.  From  (40)  we  solve  x  in. terms  of  u, 


K- 


k 


^(u2-s2) 

x  — 7 - — - 

vdsn 

Substituting  (46)  into  (45)  we  get 
2qzav 


(46) 


dsn  ^  /  dsn  y  vdsn 


•(u2-s2) 


•  u)  du  (47) 


Where  d  is  the  reduced  potential  at  the  drain,  and  is  given  by 


(Pucel,  1975) 


(48) 


Equation  (44)  for  small  <*  and  (47)  for  large  ot  are  dependent  on 
the  gate  bias  voltage,  their  transfer  characteristics  as  functions 
of  (V  a+<£)/W00  »  s2  are  shown  in  Pig.  2.  Por  the  purpose  of 

®  i 

numerical  calculation,  it  is  specified  that  <*«1 .67X10^  cm  for 
small  ot  and  of-I.SXIO4,  cm”1  for  large  c*  .  Thus  for  example,  for 
c*— 1.67X105  cm”1,  N(x)  decreases  from  lO1^  cm"5(source)  to 
9X1016  cm”5  (drain),  or  for  ««*  1.5X10*  cm"1,  N(x)  increases  from 
101?  cm"5(source)  to  1.9X10^  cm"5(drain).  The  gate  length  is 
assumed  to  be  0.6  jx.  and  the  drain-source  bias  voltage  is  assumed  to 
be  high  enough  to  ensure  saturation,  i.e.,  Vds-WQQ. 

Pig.  2  shows  significant  improvement  in  linearity  when  ot 
is  large.  However,  when  <*^0,  channel  pinches  off  quickly  at  low 
gate-source  bias  voltage.  The  reason  why  it  pinches  off  so  fast 


is  that  both  the  channel  opennin^and  the  carrier  density  are 
decreasing  functions  of  x  (see^(38)).  Therefore,  the  voltage  swing 
is  very  much  limited  and  devices  with  ®<^0  are  of  little  use. 


When  <xis  small,  there  is  no  signif leant  difference  in  the 


transfer  characteristics  when  compared  with  uniform  doping  (<x» 0) 
as  it  should  be. 


17.  DEVICE  PARAMETERS 


The  small-signal  parameters  for  the  case  of  linear  doping 
is  presented  in  this  section.  The  parameters  of  a  typical  PET  is 
assumed  to  he  z*500yu,  a®0.2 yu,  t* 0.6  Nq-101^  cm"^,  Ea*4.44  K7/cm, 

Sr»12.5,  ^1q»4500  cm2/V-sec.,  $«0.8  7,  and  va*p0Es.  (0»te*f£0,9)  • 

A.  Transconductance  g^ 

The  trans conductance  g^  is  defined  as 


v  if 


gs 


-( 


hi  ds 


si  ad 

^-av^ 


'ds 


(49) 


ds 


When  a  is  small,  the  transconductance  can  he  obtained  hy  using  (44) 
and  (49)  and  is  given  hy 


a 

(h/a)-s 


4ZTsaH0 

“W55T 


(1+*ty2) 


(50) 


When  «  is  large,  equation  (47)  can  he  numerically  differentiated 
according  to  (49).  The  transconductance  for  small  and  large  oc are 
shown  in  Pig.  3.  As  expected,  for  small  o^the  variation  of  g^  with 
gate  bias  voltage  approaches  to  that  o*  the  uniform  doping.  However, 
as  |<xi|  approaches  unity  there  is  significant  improvement  in  linearity 

°f  8b- 

B.  Gate-to-source  capacitance  Cg3 

According  to  Pucel  (1975),  the  total  charge  on  the 
gate  electrod  is  given  hy 


Where  E  is  the  y-component  of  the  electric  field,  defined  as 


E  - 

y  ay 


2W, 


2^(h/a)(1  +  o4x) 


(52) 


y**o 


Differentiating  (14)  with  respect  to  y,  Ey  hence  Qg  can  he  evaluated 

id  as  (Pucel,  1 

3Q<r  /^Qg  Ss 


cgs  is  define<i  35  (I^cel»  1975) 


«S'5V|  'Vis 


v, 


(53) 


ds 


V. 


ds 


For  small  oi ,  Q_  can  he  obtained  from  (51)  and  (52), 

O 


‘g 


2€r*0VOOz3^  ret  22) 


(54) 


And  from  (53),  Cg3  becomes 


ga 


~T 


(55) 


h/a=s 


For  large  « »  Qg  becomes 

2€r£0W00z 


‘g 


f- 


tfx)  (h/a)dx 


(56) 


Where  h/a  is  given  in  (41).  Substituting  (56)  into  (53),  Cgg  for 
large  tf  can  numerically  be  evaluated.  Fig.  4  shows  that  the  rate 
of  change  of  Qg  with  (Vg3+0^QO  is  larger  for  small  o(  than  that 
for  large  o(.  This  implies  that  C^g  decreases  as|o<j  increases,  as 
shown  in  Fig.  4. 

C.  Cut-off  frequency  F^ 


The  unit  gain  cut-off  frequency  is  given  by 


14 


T" 


®m 


(57) 


It  is  shown  in  Fig.  5  that  for  «>0  does  not  vary  with  7gs  and 
p_  f0r  large  o<  (where  xB.  »0. 9,  as  an  example)  is  almost  twice  as 

J# 

large  as  for  small  and  vanishing  <*.  This  improvement  in  can 
be  accounted  for  by  the  increase  in  g^j  and  decrease  in  Cga  due  to 
the  non-uniformity  of  doping. 
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V.  CONCLUSIONS 


Analytical  solution  of  potential  in  MESFET's  with  arbitrary 
doping  profile  have  been  presented.  The  linear  doping  profile  has 
been  treated  as  a  special  case,  in  detail.  Numerical  results  on 
device  parameters  for  linear  lateral  doping  profile  are  presented 
and  compared  .with  those  for  uniform  doping.  It  is  shown  that  there  is 
significant  improvement  in  linearity  as.  |c*/(  approaches  toward 
unity.  The  transconductance  of  the  PET  is  found  to  be  larger  for 

than  that  for  *£*0.1.  This  improvement  in  gjjj  can  be  accounted 
for  by  the  increase  in  carrier  density  in  the  channel.  Significant 
improvement  in  P^  can  be  realized  for  FET’s  with  increasing  doping 
density  from  source  to  drain(i.e.  ,<x^0) .  The  analytical  solution 
presented  here  is  sufficiently  general  and  it  can  be  applied  to 
other  types  of  profiles  such  as  exponential,  power  law  and  step 
for  either  vertical  or  lateral  doping  or  their  combinations.  In 
future  these  topics  will  be  investigated.  Although  experiments 
on  vertical  non-uniform  doping  have  appeared,  none  has  been  available 
on  lateral  or  general  doping.  Our  theory  predicts  that  the  device 
performance  depends  on  doping  profile.  Optimum  lateral  and  vertical 
doping  will  be  shown  to  be  important  for  PET  power  amplifier. 
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LIST  OP  FIGURES 


1.  Cross-sectional  diagram  of  MESF3T  showing  geometrical  dimensions 
(a=0.2^i,  Jl~ 0.6  jjl,  z=500  jx)  . 

2.  Transfer  characteristics  for  uniformly  (^=0)  and  non-uniformly 
(oc^O)  doped  GaAs  MESFET  (a=0.2  jxt  J*0.6  ji,  z» 500 yi,  Nq-1017  cm"3* 

V"2-895  V*  Vds“wOO>* 


3.  Trans  conductance  vs.,  normalized  gate  voltage  for  uniformly  (o<*  0) 
and  non-uniformly  (<*#0)  doped  GaAs  MESFET  (a=0.2  yx,  &  «0.6 jxt 
z=*500  jxt  N0»1017  cm"3,  Vda»W00«2.895  V). 


4.  Cgg  and  vs.  normalized  gate  bias  voltage  for  uniformly  (<* «0) 

and  non-uniformly  (o<j£0)  doped  GaAs  MESFET  (a*0.2  jx,  0.6 y l , 

z»500 yi,  Nq-1017  cm"3,  V^Wqq-2.895  V,*r»12.5). 

5.  Cut-off  frequency  F^  vs.  normalized  gate  bias  voltage  for 
uniformly  (o(» 0)  and  non-uniformly  (d<^»0)  doped  GaAs  MESFET 
(a=0.2  yu,  >^=0.6 jx,  z»500 yi,  ^da^OO3'^*8^3  V). 
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Telocity  Overshoot  of  Electrons  in  GaAs 
with  Space  Charge  and  Non-Uniform  Field 

By 

Ynn-chu  Wang  and  Zhzng- an  Qro1) 

Simple  analytical  transport  model  is  used  for  self-consistently  calculating  velocity  overshoot 
under  constant  current.  It  is  shown  that  space  charge  and  non-uniform  field  drastically  reduce  the 
overshooting  velocity  and  greatly  elongate  its  distance  of  effectiveness. 

Bin  einfaches  analytiachee  Transportmodell  wird  for  eine  selbstkonaistente  Bereohnung  dee  Qe- 
sehwindigkeitsaberschuaaes  bei  konstantem  Strom  benutzt.  Es  wild  geseigt,  daB  Raumladung 
and  ein  ungleichformiges  Veld  den  Geschwindigkeitsftberschufl  drastisch  redazieren  and  seine 
Wirkungsdistanx  weit  ausdehnen. 

1.  Introduction 

In  the  recent  development  of  high-speed  GaAs  logic  and  high-frequency  analog 
devices,  interests  are  focused  in  the  possibility  of  overshooting  the  electron  drift 
velocity  in  a  sufficiently  short  submicron  device.  When  space  charge  limited  current 
is  also  considered,  the  so-called  ballistic  transport  device  has  become  popular  among 
researchers.  Previous  papers  dealing  with  uniform  field  and  therefore  neglecting  the 
space  charge  effect  obtained  at  room  temperature  rather  high  overshooting  velocities 
(5  to  7.8  x  10T  em’/s)  within  distances  froln  0.2  to  0.4  pm  for  a  constant  field  ranging 
from  10  to  28  kV/cm  [1  to  4],  Since  in  small  devices,  electron  density  varies  signifi¬ 
cantly  close  to  the  cathode  contact,  space  charge  and  non-uniform  field  have  to  be 
considered  self-consistently.  Also  important  is  the  fact  that  capacitive  ourrent  in 
small  device  in  transient  state  must  be  accounted  for.  Previous  work  [8]  dealt  with 
non-uniform  field  non-self-conaistently  and  the  obtained  results  are  difficult  to  compare. 


2.  Analysis 

We  will  use  the  simple  analytical  transport  model  [4]  to  describe  the  electron  transient 
behaviour.  In  this  monoenergetic  model  approximate  energy-related  momentum  and 
energy  relaxation  times  are  assumed  to  be  valid  under  the  condition  that  the  electron 
energy  is  initially  (at  i  =  0)  maintained  at  a  level  higher  than  the  optical  phonon 
energy  Amc  in  order  to  avoid  the  unphysical  negative  energy  relaxation  time  [6].  This 
initial  level  is  found  to  be  4 hw,  =  0.148  eV,  and  at  this  level  the  energy  loss  rate 
per  unit  time  gives  at  t  =  0  the  correct  value  of  zero  before  excitation  and  correct 
values  after  excitation  for  higher  energy.  The  normalized  balance  equations  are  given 
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Yxs-chtj  Waso  and  Zmoro-asr  Qxv 


by 


gj  =  <hQ(‘)  E  -  V(«)  ® . 


(1) 


de 

dt 


a^vE  —  b^gie) , 


(2) 


where  t>,  s,  and  E  are,  respectively,  the  average  electron  drift  velocity  of  electrons, 
energy,  and  field.  The  detailed  description  of  the  energy  functions  has  appeared 
elsewhere  [4].  In  these  equations  the  (thermal)  diffusion  effects  are  assumed  to  be 
negligible  because  diffusion  is  a  slow  process  in  a  transient  problem  such  as  the  one 
treated  here. 

The  space  charge  effect  is  governed  by  the  Gauss  equation 


8  E  g 
8*  “s. 


(3) 


where  g  is  the  charge  density  and  e,  the  dielectrio  constant.  Here  we  consider  an 
undoped  GaAs  device.  The  total  current  is  given  by 


J(*)  =  «.  ^  + 


(4) 


Substituting  (3)  into  (4)  we  get 

d*  '  J(fl 

df  s,  *  (5) 

where  the  total  derivative  d/di  =  8/81  -f  ®  8/8*. 

Equations  (1),  (2),  and  (6)  constitute  the  three  governing  equations  provided  J(t) 
is  prescribed.  The  initial  conditions  are  as  follows:  v(0)  *=  E( 0)  0  and  s(0)  =  0.148. 

The  current  is  assumed  to  be  J(l)  =  J,(l  +  «  cos  cot)  6(1),  where  6(1)  is  the  unit-step 
funotion,  J,  the  bias  current;  and  *  the  ao  amplitude.  The  velocity  overshoot  results 
will  be  given  as  a  function  of  *  or  1  with  total  ourrent  J,  and  et  as  parameters.  The  ao 
frequency  of  the  excitation  current,  is  assumed  to  be  co  =  2 nf  where  /  *»  100  GHz. 
The  temperature  is  assumed  to  be  300  K. 


3.  Numerical  Results 

Fig.  1  shows  the  velocity  overshoot  for  given  bias  current  J,  versus  length  for  <x  =  0. 


Fig.  1.  Velocity  overshoot  vs.  distance  (length)  for  do 
current  excitation  only  (a  —  0).  The  bias  de  current  . 
density  is  a  parameter 
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Fig.  2.  Field  v*.  distance  for  a  «=  0  with  a*  a  parameter 


When  compared  with  the  velocity  overshooting 
corvee  for  constant  field  [4]  it  is  evident  that  the 
•pace  charge  field  causes  (i)  lower  peak  velocities  and 
that  (ii)  peaking  occurs  at  much  longer  time  and  dis¬ 
tance.  Thus  a  submicron  device  will  have  much  larger 
transit  time  than  previously  anticipated  because  the 
peaking  velocity  is  reduced  to  (2  to  4)  x  10*  om/s,  if 
the  space  charge  and  capacitive  current  are  properly 
accounted  for,  as  is  presented  here.  Collision  will  re¬ 
main  dominant  even  in  submioron  devices  due  to 


space  charge  and  non-uniform  field  [7]. 

Fig.  2  shows  the  non-uniform  field  versus  length  with  the  total  ourrent  aa  a  para¬ 
meter  for  a  mm  0.  Fig.  3  shows  the  v-E  curves  for  given  device  length  for  *  «=  0.  For 
long  devices  the  v~E  behaviors  are  very  dose  to  the  equilibrium  curve,  whereas  for 
short  devioes  they  deviate  from  the  steady  state  due  to  the  effect  that  transit  times 
for  short  length  become  comparable  to  energy  relaxation  time  [8].  There  is  a  striking 
difference  between  this  and  its  counterpart  with  constant  field  [8],  i.e.,  velocities 
decrease  drastically-beyond  the  peaks  when  a  nota-uniform  field  is  considered  as  shown 
in  Fig.  3,  but  they  become  almost  saturated  beyond  the  peaks  for  the  constant  field 
case. 


In  Fig.  4  the  total  ourrent  is  plotted  against  the  voltage  for  different  deviee  lengths 
for  *  =*  0.  For  a  0.2  pun  device  the  ourrent  behaves  like  half-power  function  of 
voltage  at  low  current,  it  increases  fast  at  medium  ourrent  levels  and  then  its  increase 
slows  at  high  ourrent  level.  For  a  1  pm  device  the  situation  is  quite  different  in  that 
the  current  increases  more  slowly  and  tends  to  saturate  at  high  ourrent  level.  Fig.  ff 
presents  the  length  versus  transit  time  relationship. 

Fig.  6  to  10  give  the  corresponding  curves  fora  *=  1.0.  In  this  case  the  sc  current 
affects  the  velocity  overshoot  drastically  as  evident  from  comparison  of  Fig.  6  and  1. 


!%•*  Fig.  4 

Fig.  3.  Velocity  va.  field  for  variable  deviee  length  and  a  —  0 
Fig.  4.  do  current  density  J%  vs.  voltage  U  for  a  «  0  and  different  deviee  length 
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Fig.  8  Fig.  6 

Fig.  5.  Device  length  v*.  transit  time  for  a  «=  0 

Fig.  8.  Velocity  overshoot  vs.  transit  time  for  both  do  and  ec  current  excitation  (a  ■=  1.0)  with 
ae  a  parameter 

Non-uniformity  of  field  versus  transit  time  as  illustrated  in  Fig.  7  shows  that  non¬ 
uniformity  of  field  is  greatly  reduced  by  the  ac  current  excitation  characterized  by  at 
and  disappears  at  low  bias  current  Jt.  Fig.  8  shows  that  ao  current  enhances  the  peak 
velocity  for  short  devices.  For  a  0.4  fim  device  the  peak  velocity  is  about  4  x  10T  cm/s 
at  a  threshold  field  of  10  kV/cm.  These  values  are  larger  than  those  for  the  case  of 
*  =  0  (Fig.  3).  Fig.  9  shows  the  bias  current  Jt  as  a  function  of  voltage  for  different 
device  lengths.  When  it  is  compared  with  Fig.  4  (for  at  =  0),  we  see  that  there  the 
voltage  drop  across  the  device  corresponding  to  the  same  bias  current  is  smaller  due 
to  the  ac  current  contribution.  In  Fig.  10  length  versus  transit  time  iB  shown. 


Velocity  Overshoot  of  Electrons  in  GaAs  with  Space  Charge 
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Fig.  9 

Fig.  9.  do  current  density  J,  vs.  voltage  for  «  «=  1.0 
Fig.  10.  Length  vs.  transit  time  for  a  «=  1.0  and  different  J% 

4.  Conclusions  and  Remarks 

In  conclusion  we  have  found  that  for  a  realistic  analysis  of  small  devices  for  which 
both  the  space  charge  and  non-uniform  field  effects  are  accounted,  the  average 
velocity  of  electrons. overshoots  only  mildly  at  rather  large  distance  when  there  is 
no  ac  current  excitation  («  =  0),  contrary  to  that  for  the  constant  field  ease.  Thus 
collision  remains  important  for  submicron  devices.  In  this  case  the  U-E  character¬ 
istics  also  show  different  behavior  in  that  though  the  velocity  increases  with  decreas¬ 
ing  length  due  to  non-equilibrium,  negative  resistance  regions  exist.  The  J9-U  curves 
show  a  great  deal  more  saturation  tendency  than  predicted  by  the  “baUiatio”  theory. 
In  the  case  of  do  plus  ac  current  excitation  (<x  «=  1.0),  the  ao  current  greatly  modifies 
the  field  to  such  an  extent  that  the  field  becomes  more  uniform  for  decreasing  bias 
current. 

It  is  beneficial,  before  this  work  is  concluded,  to  remark  on  the  special  choice  of 
0.148  eV,  instead  of  0  eV,  as  the  initial  electron  energy  at  t  =  0.  As  is  presented  in 
[4],  the  energy  loss  rate  is  dominated  by  longitudinal  polar  optical  mode.  Unfortunat¬ 
ely,  relaxation  times  valid  for  all  electron  energies  do  not  exist  for  polar  optical  modes 
due  to  partial  inelasticity  of  scattering,  although  they  are  very  useful  in  describing 
the  transport  properties  [9].  For  e-  or  T,  0D,  where  s_  =  hca,  is  the  longi¬ 
tudinal  polar  optical  phonon  energy  ati  —  0,  T,  the  equivalent  electron  temperature, 
and  0D  the  Debye  temperature  (418  K  for  GaAs),  the  scattering  is  elastic  so  that 
relaxation  times  are  well  defined.  Our  energy  relaxation  rate  formula  [4]  obtained 
from  solving  Boltzmann’s  equation  by  use  of  a  drift  Maxwellian  distribution,  whereas 
valid  for  large  energy,  does  not  hold  for  e  <  er  (near  *  =  0)  [9].  In  this  low-energy 
region  the  formula  without  correction  factor  would  give  (s  —  ejJ  r~l  ~  In  (4*/*,)//*  -» 
— •  —oo  so  that  r~l  -*  oo  as  e  -*  0,  a  result  which  so  severely  overestimates  the  contri¬ 
bution  of  the  low-energy  electrons  to  the  energy  loss  rate  as  to  cause  infinite  energy 
gain.  In  actuality  electrons  with  t  cannot  emit  phonons,  thus  they  should  not 
contribute  to  spontaneous  emission  ana  cannot  gain  energy  either.  To  circumvent  the 
unphysical  energy  gain  difficulty  we  postulate  that  the  minimum  energy  of  electrons 
in  our  simplified  model  should  at  least  be  e^,  the  thermal  lattice  energy,  i.e.,  before 
entering  from  the  oathode  into  the  GaAs  medium,  an  electron  should  have  certain 
thermal  energy  although  the  field  there  is  zero  (e.g.,  n*nn+  structure).  To  satisfy 
the  non-emission  condition  for  low-energy  electrons,  we  further  assume  that  s  >  «»  > 
>  ij,  at  300  K  and  require  that  (e  —  tj)  t~x  be  zero  at  the  cathode.  This  demands  a 
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correction  factor  such  that 

(•) 

To  fix  the  constants  r  and  »  we  follow  Ridley’s  arguments  [6],  i.e.,  they  should  be 
so  chosen  that  when  s  becomes  large,  correct  energy  loss  rate  is  obtained.  He  suggested 
that  r  =  4  and  n  =  2.  The  initial  energy  is  therefore  taken  to  be  4«p  =  0.148  eV 
instead  of  the  commonly  assumed  0  eV  in  typical  Monte-Carlo  calculations. 

A  possible  physical  explanation  is  as  follows:  Although  initially  the  electron  has 
no  drift  velocity,  v(0)  —  0  since  E(0)  =  0,  and  therefore  no  drift  energy,  continuity 
of  current  however  requires  the  diffusion  current  to  be  present  at  the  cathode  where 
an  electron  temperature  gradient  exists;  the  cathode  does  not  require  the  presence 
of  infinitely  large  electron  density  as  commonly  assumed.  This  implies  that  a  thermal 
spread  of  the  cathode  exists,  inside  of  which  the  electron  temperature  must  be  higher 
than  300  K.  As  soon  as  the  electron  enters  the  medium,  it  picks  up  the  drift  energy 
at  the  high  rate  proportional  to  »*;  thus  the  further  it  moves  away  from  the  cathode, 
the  more  predominant  the  drift  energy  is,  until  it  reaches  the  anode  (which  we  do 
not  consider  here),  where  the  velocity  and  field  again  vanish  and  the  diffusion  current 
takes  over  as  the  total  particle  current  as  before.  It  may  be  added  that  the  diffusion 
current  arises  from  the  fact  that  an  electron,  shortly  before  entering  the  anode,  has 
very  high  temperature  gained  from  the  field,  thus  a  temperature  gradient  exists  at 
the  anode. 

Once  the  initial  energy  is  properly  determined,  the  detailed  diffusion  process  on 
the  boundary,  where  electrons  leave  or  enter,  can  be  disregarded  in  our  simplified 
model.  It  remains  to  be  explained  as  to  why  diffusion  can  be  neglected  in  the  n-GaAs 
region.  It  is  not  difficult  to  include  diffusion  terms  in  the  balance  equations  (1)  and 
(2),  but  the  transient  analysis  becomes  tedious.  Instead  of  doing  thiB,  one  can  include 
the  diffusion  current  JD  in  the  phenomenologic  current  equation,  which  is  somewhat 
easier  in  getting  order-of-magnitude  comparison.  JD  is  proportional  to  the  gradient 
of  density  q  governed  by  (3)  and  the  field-dependent  diffusion  coefficient  D(E).  To 
a  firat  approximation.  Fig.  2  for  dc  current  excitation  only  shows  that  the  field  E  is 
almost  linear;  thus  8g/ 8x  is  not  large,  justifying  negligibility  of  7S.  Fig.  7  shows  that 
for  dc  plus  ac  excitation  q  does  indeed  vary,  but  for  the  high-field  region  D(E)  decreases 
drastically  and  JB  still  may  be  small  enough  to  be  negligible.  It  is  noticed,  however, 
that  near  the  cathode,  where  E  is  small  and  therefore  3(E)  large,  diffusion  current 
may  generally  become  important.  But  even  for  this  situation  our  special  initial  con¬ 
dition  seems  to  work  to  our  advantage  giving  small  gradient  of  density  thus  favoring 
drift  current. 
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As  reported  in  the  interm  reports  we  have  had  several  cou¬ 
pling  activities  with  other  government  laboratories.  With  Harry 
Diamond  Laboratories  we  have  fabricated  ohmic  contacts  on  super¬ 
lattice  structures.  With  the  Naval  Reasearch  Laboratories  we 
have  been  investigating  the  annealing  and  characterisation  of 
high  energy  implants.  In  addition,  we  have  been  investigating 
DLTS  spectra  of  irradiated  materials  .  One  of  these  irradiation 
studies  has  been  concluded  and  a  paper  is  attached. 
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b  Abstract  i 

•  •  *,  ..  —  A 

Deep  Level  Transient  Spectroscopy  (DLTS)  has  been  used  to  study  the  , 

nature  of  the  neutron-irradiation  induced  defects  In  GaAs.  The 
measurements  were  made  on  a  group  of  Schottky  barrier  samples  on  n-type 
(1,2X10*'  cm*3)  LPE  GaAs.  The  samples  were  exposed  in  a  reactor  to  •  •  . 
neutrpn  fluences  from  5X1014  cm”2  to  5X1016  cm”2.  DLTS  measurements 
on  unirradiated  samples  revealed  no  electron  traps  with  concentrations 
greater  than  2X1013  cm"3.  The  samples  irradiated  to  5X101*  and  _• 

2X1015  cm"2  were  studied  as  Irradiated  and  after  annealing  to  temperatures 
of  500C.  The  positions  of  the  peaks  studied  here  are  sensitive  to  the 
the  electric  field  in  the  depletion  layer.  Evidence  Is  also  * 
presented  to  t'amonstrate  that  the  strength  of  the  electric  field 

enhanced  emission  effect  decreases  on  annealing.  The  spectra  V  , -••-v 

for  samples  annealed  to  500C  contained  three  lines  at  .38, .56  and  .78  eV  J»  .  _ 

below  the  conduction  band.  The- line  at  .78  eV  was  the  most  prominent 
of  these,  and  Its  line  shape  was  found  to  within  about  10X  of  that  expected 
for  a  point  defect. 


\ 
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Introduction 


Following  the  observatlon.by  Electron  Paramagnetic  Resonance  (EPR)  that 
the  Asga-  antis, Ite  In  GaAs  could  be  Introduced  by  neutron  Irradiation  (1)  * 
there  were  several  attempts  to  determine  whether  the  Deep  Level  Transient 
Spectroscopy  (DLTS)  line  EL2  is  also  produced  by  neutron  Irradiation  (2,3,4). 
G.  M.  Martin  and  S.  Makrara-Ebeld  reported  observing  that  neutron-irradiation 
produced  a  broad  band  and  some  of  Its  changes  on  annealing  (2,4). 

The  general  conclusion  from  these  experiments  was  that  the  Introduction 
rate  for  EL2  was  much  smaller  than  the  Introduction  rate  for  the  ASg- 
antl site ‘found  by  EPR.  The  main  feature  of  the  DLTS  spectra 
was  a  broad  line  centered  about  -40C  with  some  small  shoulders.  A  broad 
line  like  this  Is  assumed  to  be  associated  with  complex  defects  (5)  rather  j 
than  a  well  defined  point  defect.  .  i 
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There  are  a  number  of  reasons  for  making  the  In-depth  DLTS  study  of  the 
broad  neutron-irradiation  Induced  line  which  Is  reported  here.  One  Is  the 
fact  that  the  Introduction  rate  (l)*for  this  line  is  near  to  the  rate  found 
In  the  EPR  experiments  on  the  Asc,  antisite.  This  becomes  even  more  . 
compelling  when  It  is  understood  that  the  nature  of  the  electron 
wavefunctlon  associated  with  the  EPR  signal  Implies  that  the  EPR  experiments 
are  not  very  sensitive  to  whether  the  antisite  Is  complexed  with  another 
defect  or  a  donor.  This  broad  line  Is  also  the  dominant  electron  trapping* 
feature  observed  In  the  upper  half  of  the  band  gap  In  DLTS  experiments. 

Liquid  Phase  Epitaxial  (LPE)  samples  were  chosen  for  this  experiment 
because  It  Is  possible  to  grow  LPE  GaAs  which  is  essentially  free  of 
electron  traps.  This  eliminates  the  need  to  subtract  background  signals 
from  defects  present  In  the  starting  material.  The  nature  of  the  radiation 
Induced  defects  may  also  be  changed  by  Interactions  with  defects  already 
present  In  the  starting  material.  Heavily  doped  (1.2X1017  Sn  cm"3)  layers 
were  chosen  In  order  to  use  a  range  of  fluences  similar  to  those  used  In  the 
EPR-  studies.  At  least  1X1015  neutrons  cm'2  are  needed  to  obtain  a 
reasonable  signal  to  noise  ratio  In  EPR  work,  while  fluences  much  higher 
than  this  are  expected  to  compensate  even  heavily  doped  GaAs  making 
.conventional  DLTS  experiments  difficult. 

Experiment 

Th# samples  used  In  this  experiment  were  250  micrometer  diameter  gold 
Schottky  diodes  formed  on  a  10  micrometer  thick  n-type  (1.2X1017  sn  cm-3) 

LPE  grown  GaAs  layer.  Several  GaAs  chips  were  prepared  each  with  about  ten 
diodes.  Each  chip  was  Irradiated  In- a  reactor  with  a  different  neutron 
fluence  to  obtain  a  group  of  samples  with  fluences  between  5X10l*  an^ 

5X101®  neutrons  cm"?.  The  neutron  Irradiations  were  carried  out  by  E. Weber 
In  a  manner  similar ‘to  that  used  by  him  to  prepare  samples  for  EPR 
experiments.  The  thermal  annealing  treatments  were  carried  out 
In  a  hydrogen  atmosphere  after  stripping  off  the  gold  Schottky  metallzatlon. 
Only  those  samples  with  fluences  of  5X1014  and  2X1015  cm"2  are  reported  here 
as  the  others  were  rendered  semi -Insulating  by  the  Irradiation. 


.  The  DLTS  measurements  were  carried  out  with  a  one  MHz  capacitance' 
meter  by  applying  a  quiescent  reverse  bias  and  then  periodically  applying 
a  voltage  pulse  to  fill  the  electron  traps.  A  computer  controlled  data 
acquisition  system  was  used  to  set  the  temperature  and  to  digitally  record 
the  capacitance  transient  after  each  filling  pulse.  The  difference  between 
the  capacitance  at  times  Tj  and  T2  was  determined  by  a  computer  program  to 
simulate  a  boxcar  spectrum.  The  spectra  generally  were  recorded  from  liquid 
nitrogen  temperatures  to  about  420  K.  C-V  data  were  measured  at  several 
temperatures  to  obtain  the  capacitance  at  the  quiescent  reverse  bias  used 
during  the  DLTS  measurements  and  to  calculate  the  free  carrier  concentrations 


Results 


The  quality  of  the  LPE  layers  was  checked  prior  to  the  neutron 
Irradiations  by  measuring  DLTS  and  C-V  data  which  Indicated  that  the 
•  electron  trap  concentration  was  loss  than  2X1013  cm*3.  Layers  with  carrier 
concentrations  of  1.2X1017  cm"3  were  chosen  in  order  to. use  neutron  . 
fluences  near  those  used  In-  EPR  experiments  on  the  As^  antisite. 


Simulated  boxcar  spectra  for  samples  irradiated  to  fluences  of  5X10^ 
and  2X101&  neutrons  cm"2  are  Illustrated  in  Figs'.  1  and  2  respectively 
with  Ti»5ms  and  T2*50ms.  The  OLTS  spectra  shown  in  Fig.  la  Is  typical  of 
those  often  observed  for  neutron  Irradiated  GaAs  In  that  It  Is  very  broad. 


to  zero  bias;  c:  bias  «  3  volt  pulse  height  1  volt. 


The  full  width  at  half  the  peak  amplitude  for  this  line  Is  about  150K  which 
Is  four  to  five  times  broader  than  the  lines  observed  for  point  defects.  The 
width  Is  a  sign  that  a  complex  defect  Is  Involved  In  the  trapping.  These  data 
were  obtained  with  a  3  volt  reverse  bias  and*  by  pulsing  to  zero  bias  to  fill 
the  traps.  As  a  furthur  study  of  the  sample  the  spectra  In  Fig.  lb  was 
recorded  with  a  one  volt  bias  and  a 'zero  bias  filling  pulse.  The  peak  In 
Fig.  lb  Is  at  a  higher  temperature  than- that  In  Fig.  la,  suggesting  either 
a  distribution  of  defects  which  varies  with  depth  1.n  the  sample  or  that  the 
emission  rate  In  Fig.  la  was  enhanced  by  the  larger  average  electric  field 
used  In  that  measurement.  To  determine  which  effect  Is  responsible  for  the 
peak  shift,  the  data  In  figure  lc  were  measured  with  a  3  volt  reverse  bias 
and  a  one  volt  high  filling  pulse  to  obtain  the  same  average  electric  field 
used  In  Fig.  lb,  but  over  a  different  sample  volumne.  Since  almost  Identical 
spectra  are  found  In  Figs,  lb  and  lc  the  conclusion  Is  that  the  large  average 
electric  field  used'for  the  data  In  Fig.  la  enhanced  the  emission  rate  and 
caused  the  peak  to  shift  to  lower  temperatures. 

An  estimate  of  the  number  of  defects,  N*,  can  be  made  by  using  I 


Nt»FAC„  Ns/Cp  .  (1)  .  : 

I 

ACp  Is  the  peak  amplitude,  Cq  the  capacitance  at  the  quiescent  reverse 
bibs,  Ns  Is  the  free  carrier  concentration  and  F  Is  a  constant.  Poisson's  ! 
equation  can  be  used  to  evaluate  F  for  the  case  of  a  point  defect  with  an 
exponential  decay.  In  general  F  depends  on  the  bias,  pulse  helght.bullt  In 


potential,  defect  energy  as  well  as  T}  and  T£.  For  the  various  conditions 
used  In  these  experiments  the  approplate  values  of  F  ranged  from  5  to  20. 

Using  these  assumptions  results  In  a  calculated  defect  concentration  of 

2X10*5  cm”3  for  Fig.  1  which  corresponds  to  a  defect. introduction  rate  of  4. 

• 

This  should  be  an  underestimate  of  the  actual  value  because  the  analysis 
Ignores  the  fact  that  the  line  results  from  the  distribution  of  defects. 

Since  the  entire  capacitance  transient  was  recorded.  It  Is  possible  to  measure 
Ct,  the  amplitude  directly  after  the  filling  pulse,  and  to  use  It  In  place 
of ACp  In  Equation  1  along  with  an  appropriate  value  of  F  to  find 
Nt"7X1015  cm"3.  This  gives  an  Introduction  rate  of  14,  which  Is  probably 
a  better  estimate  of  the  Introduction  rate  as  It  represents  an  Integration  of 
all  the  transients  under  the  peak. 


Data  on  the  annealing  behavior  of  a  sample  Irradiated  to  a  fluence 
of  2X1015  neutrons  cm"2  are  presented  In  Fig.  2.  The  unannealed  sample 
Is  shown  In  Fig.  2a  along  with  the  results  of  annealing  It  at  400C,  450C 
and  500C  In  Figs. 2b, 2c  and  2d.  While  all  these  curves  were  recorded  with 
the  same  bias  and. pulse  height,  measurements  were  also  made  with  a  number  of 
other  bias  and  pulse  height  combinations.  The  results  of  those  tests  Indicated 


that  the  emission  rates  were  electric  field  enhanced  for  all  but  the  sample 

annealed  at  500C.  The  data  shown  In  Fig.  2  are  evidence  for  a  weakening 

of  the  processes  responsible  for  the  electric  field  enhancement  effects. 

The  electric  field  applied  during  these  experiments  Increased  on  annealing 
because  the  depletion  layer  width  decreased  due  the  free  carrier  concentration 

recovering  from  3.6X101®  to  9.5X101®  cm“^.  The  arrow  on  each  of  the  curves  In 

Fig.  2  Indicates  the  position  of  the  peak  of  the  spectra  and  they  show  that 
the  peak  Is  slowly  moving  to  higher  temperatures  on  annealing.  This 
signifies  a  decrease  In  the  emission  rate  at  a  given  temperature; 

A  furthur  examination  of  Fig.  2  Indicates  the  growth  of  a  shoulder  on 
the  low  temperature  side  of  the  400C  and  450C  data.  A  high  temperature 
shoulder  was  also  found  on  the  4S0C  data  though  It  cannot  be  seen  In  Fig.  2. 
The  low  temperature  shoulder  did  not  exhibit  any  evldencfe  of  electric  field 
dependence.  The  high  temperature  shoulder  was  not  distinct  enough  to 
determine  If  It  was  field  dependent. 

The  most  significant  change  In  the  spectra  occured  after  annealing 
at  500C.  The  Intensity  of  the  lines  decreased  significantly  and 
a  new  line  appeared  at  the  high  temperature  end  of  the  spectra.  This 
was  the  first  observation  of  a  peak  this  high  in  temperature  In  these 
samples.  The  shape  of  the  line  agrees  well  with  that  predicted  by  a  point 
defect  model.  A  similar  line  Is  also  observed  In  samples  irradiated  at 
hlgher^fluences  and  annealed  at  500C  to  restore  their  conductivity.  The 
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Fig.  3.  Arrhenius  plot  for  the  spectra  In  Fig.  2c  with 

*a  -.38  eV.  A  -3.3xl04  sec ^  JT2s  Ea  -.56  eV, 

A  -1.4x10®  sec”1  JT2:  Em  -.78  eV,  A  -1.3x10*  sec”1  JT2 
and  with  literature  values  of  E.  -.825,  A  -2.9xl07  sec” 
for  EL2  (6). 


annealing  process  was  checked  on  unirradiated  samples  to  confirm  that  the 
new  line  was  not  Introduced  during  the  processing. 

An  Arrhenius  plot  for  the  lines  In  Fig.  2c  Is  shown  In  Fig.  3  along 
with  a  line  at  the  position  expected  for  EL2.  The  data  in  this 
experiment  were  fit  to 

€  -  AT2EXP(-Ea/KbT).  (2) 


The  high  temperature  line  In  Fig.  3  has  Ea*.78  eV  and  A*1.3X107  which  are 
near  the  values  of  Ea*.825  eV  and  A*2.9X10'  often  quoted  for  EL2  (61. 

The  Intensity  of  the  .78  eV  line  corresponds  $o  a  defect  concentration 


The  Intensity  of  the  .78  eV  line  corresponds  to  a  defect  concentration 
of  1.5X10*5  cnj-3.  ^is  needs  to  be  compared  with  a  defect  concentration 
estimated  to  be  between  1.5X1016  and  2.5X1016  cm-3  for  the  unannealed  data 
In  Fig.  2a. 

Olscusslon 

The  results  presented  here  give  some  new  Insite  Into  the  nature  of  the 
broad  line  observed  In  DITS  experiments  on  neutron  Irradiated  GaAs.  The 
first  point  to  note  Is  that  the  position  of  the  peak  for  both  samples 
presented  here  Is  very  sensitive  to  the  strength  of  the  electric  field 
In  the  junction.  To  estimate  the  position  of  the  peak  at  zero  field,  the 
temperature  of  the  peak  for  the  sample  shown  in  Fig.  1  was  plotted  for 
several  emission  rates  against  the  applied  field  which  ranged  from 
1.2X10®  to  1.8X10®  volts/cm.  On  extrapolating  the  data  to  zero  electric 
field  the  lines  fell  In  the  temperature  range  where  EL2  Is  found.  While 
this  suggests  that  the  line  may  be  a  field  shifted  EL2  signal  this  Is  not 
likely  because  EL2  has  not  been  observed  to  be  sensitive  to  fields  as 
low  as  1.2X10®  volts/cm  (71.  In  addition  when  the  line  width  Is 
extrapolatd  to  zero  field  It  Is  too  large  for  EL2. 


too  large  for  EL2. 


The  line  width  suggests  that  the  line  Is  associated  with  a  complex 
defect  with  a  distribution  of  energies.  To  obtain  an  estimate  of  the 
width  of  the  defect  distribution  a  program  was  developed  to  find  the 
line  shape  that  results  when  several  point  defects  with  varrylng  energies 
and  concentrations  are  added  together.  The  energies  and  concentrations  of 
the  model  lines  were  adjusted  to  obtain  a  good  fit  to  a  spectrum  measured 
at  low  reverse  bias  for  the  sample  In  Fig.  2.  A  distribution  extending 
from  about  .46  eV  to  .72  eV  with  a  full  width  at  half  amplitude  of  .1  eV 
was  found  to  give  good  fits  to  the  line  shapes  measured  for  a  variety 
of  emission  rates.  The  peak  of  the  distribution  was  near  the  energy 
obtained  In  an  Arrhenius  plot  of  the  DLTS  line. 

Another  Important  feature  of  the-  data  Is  the  decrease  in  the  strength 
of  the  field  effect  on  annealing.  This  suggests  that  the  strength  of  the 
mechanism  Is  tied  to  the  complexes  which  are  being  annealed.  The  data 
presented  here  Indicate  a  considerable  change  takes  place  on  annealing 
at  500C  compared  to  the  changes  produced  on  anneal Ingh  at  400C  and  4500. 
Because  of  this  rapid  change  between  450C  and  500C  It  was  not  possible 
to  track  the  development  of  .  the  peak  at  .78  eV  In  the  500C  data.  This 
line  may  result  from  a  weakening  of  the  field  effect  or  It  may  Indicate 
the  development  of  a  point  defect  as  the  complexes  break  up.  Since 
the  OUTS  data  still  show  a  broad  peak  at  low  temperatures  and  the  free  • 
carrier  concentration  Is  not  fully  recovered  at  5000,  there  are  a 
number  of  reasons  to  expdct  that  the  .78  eV  line  may  be  due  to  the 
defect  associated  with  EL2  complexed  with  another  defect. 


1 

.1 


Several  comparisons  may  be  made  between  these  results  and  those  ob> 
talned  by  EPR  measurements  of  the  AsGa  antisite.  The  Introduction  rate 
found  here  Is  between  4  and  12  cm-1  which  compares  well  with  the  rate  of 
10  carl  found  for  EPR  data  on  the  Asc.  anti site  In  semi -insulating 
I  samples  Irradiated  at  the  same  time  (8).  A* review  of  the  EPR  data  for 
1  the  annealing  behavior  of  AsQa  antisites  found  In  as  grown  material  and 
for  anti sites  produced  by  either  neutron  Irradiation  or  plastic  deforma¬ 
tion  suggests  that  antisites  exist  In  a  variety  of  complexes  with 
different  thermal  stabilities  (9).  The  DLTS  spectra  presented  here  show 
neutron-irradiation  Induced  defects  which  anneal  out  between  400  and 
|  500c.  A  similar  annealing  stage  has  been  observed  by  EPR  for  neutron 

irradiated  GaAs  il). 


In  summary,  the  field  dependence  of  the  broad  neutron  irradiation 
induced  DLTS  line  and  the  appearance  on  annealing  of  a  line  near  the  EL2 
.  position  suggest  that  the  Irradiation  has  produced  the  defect  responsible 
I  for  the  EL2  line  but  complexed  with  other  defects.  These  other  defects 
result  In  a  line  which  is  very  sensitive  to  .the  applied  electric  field. 

In  addition,  we  conclude  that  these  DLTS  results  have  many  of  the 
same  properties  as  the  Asga  antisite  measured  by  EPR. 
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